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Abstract 

We show how to perform universal Hamiltonian and adiabatic computing using a time- 
independent Hamiltonian on a 2D grid describing a system of hopping particles which string 
together and interact to perform the computation. In this construction, the movement of one 
particle is controlled by the presence or absence of other particles, an effective quantum field ef¬ 
fect transistor that allows the construction of controlled-NOT and controlled-rotation gates. The 
construction translates into a model for universal quantum computation with time-independent 2- 
qubit ZZ and XX+YY interactions on an (almost) planar grid. The effective Hamiltonian is arrived 
at by a single use of first-order perturbation theory avoiding the use of perturbation gadgets. The 
dynamics and spectral properties of the effective Hamiltonian can be fully determined as it corre¬ 
sponds to a particular realization of a mapping between a quantum circuit and a Hamiltonian called 
the space-time circuit-to-Hamiltonian construction. Because of the simple interactions required, 
and because no higher-order perturbation gadgets are employed, our construction is potentially 
realizable using superconducting or other solid-state qubits. 
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I. INTRODUCTION 

The first proposals for quantum computers used time-dependent Hamiltonians to enact 
unitary quantum logic gates m and the first prototype quantum computers were realized 
using such time-dependent methods via electromagnetic resonance @E|. In 1986, Feynman 
proposed a method for performing quantum computation using a time-independent Hamil¬ 
tonian [8]. His motivation was to make a model of quantum computation that resembled 
more closely the time-independent Hamiltonian dynamics of the fundamental laws of physics. 
Feynman’s trick was to adjoin a global clock variable that regulated the pace of the compu¬ 
tation. In 1987, Margolus constructed a model for Hamiltonian quantum computation that 
was spatially homogeneous, eliminating the pointer variable by embedding the computation 
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in an asynchronous cellular automaton [9], In the Margolus method, the role of the pointer 
variable is subsumed in the positions of clock particles or ‘chronons’ that carry quantum 
bits with them as they progress through the computation. The global clock variable is thus 
replaced by a local clock variable, one for each degree of freedom in the computation. 

The Feynman and Margolus models were originally conceived as dynamic models in which 
the clock or the chronons are prepared in a traveling wave state that propagates through the 
computation. Later analyses of these models have shown that it is not necessary to prepare 
the initial state in traveling wave state: one can just initialize the input state at fixed initial 
clock time (see e.g. [ID]). 

It was also noted that the Feynman Hamiltonian allows quantum computation to be em¬ 
bedded in the ground state of the quantum system; it provides a means for mapping the 
time-dynamics of a quantum system or a quantum circuit onto the ground-state of a master 
‘circuit’ Hamiltonian which includes the dynamics of the clock variable m- This ground 
state can then be reached via adiabatic quantum computation m The practical disad¬ 
vantage of the universal circuit Hamiltonian obtained through this construction is that it 
involves many-qubit interactions or qudit degrees of freedom. Such interactions can rewrit¬ 
ten in terms of simpler, say, 2-qubit interactions through the use of perturbation theory 
(by using so-called perturbation gadgets [TDffTTj h However, effective interactions obtained 
in fcth-order degenerate perturbation theory with perturbative coupling g and gap A of 
the unperturbed Hamiltonian scale in strength as g(g/ A) fc_1 leading to a correspondingly 
small gap of the effective Hamiltonian (as compared to the physical device temperature). 
In addition, multiple uses of (higher-order) perturbation theory lead to Hamiltonians with 
undesirable qubit overhead and complexity. Consequently, existing models of Hamiltonian 
quantum computation based on pairwise qubit interactions are not particularly suitable for 
physical implementation using, e.g., solid-state quantum information processors. 

The Margolus asynchronous cellular automaton model of Hamiltonian quantum com¬ 
putation relies on spatially homogeneous interactions which allow the chronons that carry 
the computation to progress at different rates at different points: this construction can be 
thought of as a quantum-computation based model of Wheeler’s ‘many-fingered time ’ ra. 
and has been recently formalized in [19] under the name space-time circuit-to-Hamiltonian 
construction. Several works have either explicitly or implicitly formulated proposals of do¬ 
ing universal Hamiltonian computation |20| 21j or quantum adiabatic computation [221125] 
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using this construction. A related proposal which seeks to do adiabatic computation using 
the idea of quantum adiabatic transistors has been formulated in [[SB], [27]. Ref. [28] has 
proposed a way of doing Hamiltonian computing using the Feynman construction using 
only 2-qubit interaction and no application of perturbation theory. However, in order to 
make this model geometrically local on a 2D grid, the author estimates that each qubit is 
involved in a high number, at most 28, qubit interactions. In addition, the use of the global 
Feynman clock leads to an inefficient usage of space (i.e., number of qubits) and time re¬ 
sources in Hamiltonian computing: the gates in the original circuit are executed sequentially 
in the Hamiltonian computation so that the time-duration of the computation scales as a 
polynomial in the size of the original circuit. In contrast, in the space-time construction the 
duration of the Hamiltonian or adiabatic computation scales with the depth of the original 
quantum circuit that the computation simulates. 

Our results are an improvement over the previous constructions in the following sense: 

The 2D grid Hamiltonian proposed in m for running an autonomous programmable 
quantum computation is based on a strong attractive interaction between hopping spin- 
1/2 particles so that in first-order perturbation theory these particles hop together through 
spatial areas where their internal states are changed according to the circuit to be im¬ 
plemented. This particular realization uses, similarly to the current work, an attractive 
interaction between particles. If one translates this construction to qubits the author’s use 
of holonomic computation to implement logic results in 4-qubit interactions. If the con¬ 
struction is described with particles with spin, it uses attractive particle interactions, single 
particle hopping terms and terms which couples the spins of pairs of particles. In addition, 
our construction goes beyond the results in Refs. [221 28] by showing how to do adiabatic 
computation using similar 2-qubit interactions and known bounded gap above the ground 
state. 

In [2EJ the resulting Hamiltonian requires particular 4-qubit interactions and thus a fur¬ 
ther use of perturbation theory. Even though we reduce particular 4-qubit interactions to 
2-qubit interactions by means of a single use of 1st order perturbation theory, it is unlikely 
that this trick is possible for arbitrary 4-qubit interactions. The Hamiltonian and adiabatic 
computation that we propose bears the closest resemblance to the 2D grid model Hamilto¬ 
nian in [25]: we show how mathematical results concerning this construction directly carry 
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over to the results in this paper. The particle interactions in [25] (and the original |24j ) 
require pairs of particles to jointly move or hop while updating their internal states : we ex¬ 
pect that a reduction of such interactions to 2-local interactions is inefficient, thus strongly 
favoring the novel construction in this paper. Two crucial steps give us this improvement 
over ea. namely the use of the railroad switch idea 0 IH2J EH] to do a classical CNOT or 
Toffoli gate, and the use of 1st order perturbation theory to coordinate the joint movement 
of the particles so that they execute the computation jointly. 

A goal of this paper is to construct models for Hamiltonian and adiabatic quantum 
computation that are compatible with experimentally available couplings in a physical device 
such 2D arrays of spin qubits in semi-conducting quantum dots or superconducting qubits. 
The recent demonstration of large-scale quantum annealers using superconducting circuits 
suggests that adiabatic or Hamiltonian computation might provide a fruitful method for 
performing quantum computation. While existing quantum annealers realize a transverse 
Ising model system, the more general couplings afforded by superconducting and other solid- 
state and quasi-solid state systems such as optical lattices might be used to perform general 
adiabatic quantum computation. We discuss possible physical realizations in Section V A| 

Running a solid-state (quantum) computation using a time-independent Hamiltonian has 
the advantage of requiring no AC fields on chip, thus removing the challenge of placing active 
control lines between quantum degrees of freedom and turning interactions on and off. A 
disadvantage of running the quantum degrees of freedom through stationary, in-place gates 
is that we trade a time-dependence for a space-dependence and thus a ID quantum circuit 
is executed on a 2D grid. 


A. Quantum Field Effect Transistor and Perturbation Theory 

An essential element in our construction is that we use high energy ‘blocking’ terms 
(Hq) to prevent lower energy hopping or kinetic terms (V) from operating. Blocking can be 
thought of as a coherent quantum version of the field-effect transistor that underlies classical 
computation: the presence of particles (electrons/holes) in the gate electrode prevents other 
particles from moving from source to drain. In classical computation the high-energy terms 
correspond to a tuning of the classical electrostatic potential by the gate electrode so that a 
sufficiently high barrier, larger than the kinetic energy of the particle at the source, blocks 
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the forward motion. In the quantum version the potential itself is formed by the presence 
or absence of a (charged) particle. Considering that the hopping particle can always tun¬ 
nel through the energy barrier, one has to choose the kinetic energy perturbatively weak 
compared to the height of the barrier. 

We illustrate the idea of blocking with a small example, shown in Fig. [l|a). Take the state 
|1) to represent the presence of a ‘particle’ on a site (represented in Fig. 1 by a filled circle 
•) and the state | 0 ) to represent the absence of a particle (represented by an empty circle o). 
Consider a Hamitonian hf 12 3 = H 0 + gV where H 0 = A|ll)(ll | 13 and V = — (cr^cr^ + cr 3 ° 2 ~)- 
If A g the term H 0 ‘blocks’ a particle from making the transition from site 2 to site 3 if 
another particle already occupies site 1 . H 123 has a simple block-diagonal form in the basis 
| 0 )( 011 and 11 )( 111 , that is 


h/123 — 


(gv 0 

\ 0 A11)(11 3 + gV 


( 1 ) 


In the 10)(01 1 part of the low-energy sector the Hamiltonian is simply gV. I 11 the |l)(l|i 
sector the eigenstates on sites 2 and 3 couple the unperturbed low-energy state 1 10)23 an d 
high-energy state |01) 23 , he \i>±) oc ^|01) 23 + (1 ± a/ 1 + 4^ 2 /A 2 )110) 23 with energy E± = 
y(l =F y/1 + Ag 2 / A 2 ). To zero’th order in g/A the low-energy state equals \i/ + ) ~ 110 )23 
while in first-order in g/A \i/>+) has an amplitude 0(g/ A) for the high-energy state |01) 2 3i 
the maximum probability of finding the system in the high-energy sector over time is thus 
suppressed to 0(g 2 / A 2 ). The probability of finding the system in the high-energy sector 
could be further decreased by interactions with the environment at temperature T < E. 
We can thus say that to lowest order in g , there exists an effective low-energy Hamiltonian 
whose dynamics is decoupled from the unperturbed high-energy sector and which equals 

H eff = slOXOli ® V - ^|l)(l|i ® |10)(10| 23 . (2) 


We see from this simple example that in the presence of blocking, 2-local hopping terms 
can take on a conditional nature: a particle can only hop from one site to another if it is 
not blocked. As a consequence, the effective low-energy Hamiltonian H e & is now 3-local. 
When performed coherently using multiple particles, by the method described below, such 
blocking allows one to build up a quantum computation out of a sequence of blocking and 
hopping moves. We note that the idea of creating blocking energy barriers for certain particle 
occupations is also used in various first-order perturbation gadgets in [29lj . 
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In our detailed construction in Section |TT] the hopping of particles is controlled by the 
presence of particles nearby. In particular, the presence of a neighbor particle provides a 
2-site energy well for another particle. This means that the energy barrier for a particle to 
go off by itself is high on all sites where there are no particles nearby and thus tunneling 
through this long energy barrier is completely suppressed. 

When the blocking site separates a ‘source’ and ‘drain’ site, as in Fig. (lb), the resulting 
device makes up a Hamiltonian quantum held effect transistor (qFET), where the particle 
at site 2 can move to site 4 only if site 3 is not blocked. Just as in a conventional electronic 
computer, where held effect transistors are used to guide charged particles through the 
computation and to perform logic, the quantum held effect transistor will be used to guide 
qubits through the computation and to perform quantum logic. The detailed construction 
of the use of quantum held effect transistors to perform Hamiltonian quantum computation 
will be given in Section V. 

Note that even when the particle at site 2 is blocked by an energy barrier at 3, it can 
still tunnel through this barrier to site 4, leading to incorrect dynamics at higher order in 
perturbation theory. To suppress such tunneling, we can add a sequence of intermediate 
sites over which the energy barrier remains high: the details of this construction are given 
in Section IIB. 

For a general analysis of the effect of blocking and weak particle hopping for many par¬ 
ticles we can invoke the results of degenerate perturbation theory. Let us briehy summarize 
the results of the systematic Schrieffer-Wolff peturbative method US [30]. Let H 0 be a 
many-particle Hamiltonian with an energy gap A between a degenerate ground-state sector 
with, say, zero energy, and higher-energy states. For example, the Hamiltonian Hq will 
add energy penalties when particles do not stay together, i.e. it is blocking a particle from 
moving away while other particles stay behind. Let Pq be the projector onto the degenerate 
ground state sector of this Hamiltonian. Let a perturbation gV be added to Ho where V 
is a sum of hopping terms for the particles and g <C A. One can show [3Uj| that there 
always exist a unitary transformation U such that U(H 0 + V)W is block-diagonal in the 
original low- and high-energy eigenspaces of Ho . One defines a ‘low-energy Hamiltonian’ 
H = P 0 U(H 0 + V)JJlPo which has then, by definition, the same spectrum as part of H 0 + V. 
One can define a perturbative expansion for the unitary operator U which gives rise to a 
Taylor expansion for H in g/A. Cutting off this Taylor expansion for H at a given desired, 
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say fcth, order in g/A then gives us an approximate effective Hamiltonian H e ^. One expects 
that this approximate effective Hamiltonian gives a proper description of the low-energy 
dynamics of U(H 0 + V)W up to errors in energy eigenvalues which scale in strength as 
g(g/A) k (and scaling extensively with system size). The unitary transformation U can be 
thought of as a dressing of the unperturbed eigenstates by the perturbation: if we develop 
the expansion up to kth. order, these dressed states are correct up to amplitudes of strength 
(g/A) k . 

To implement the blocking ideas of the quantum held effect transistor we only use per¬ 
turbation theory to lowest order so that the effective Hamiltonian equals H e " = gPoVPo, 
i.e. of strength g, and energy corrections are 0(g 2 / A). At this order U = I and the dressed 
states are the unperturbed eigenstates of Hq. In the simple example of Fig. [l] this first-order 
effective Hamiltonian is only the first term in Eq. (J2|) while in second-order one obtains the 
entire expression in Eq. (J2]). 

In our analysis we will also discuss the effect of higher-order terms in a perturbative 
expansion of H. We can show that such higher-order terms do, assuming an otherwise 
error-free running of the computation, not affect the logic of the quantum computation that 
is executed, although it will affect the spectral properties of the Hamiltonian, the form of 
the dressed states and thus the time-dynamics. 

It is important to note that errors in the perturbatively derived effective dynamics always 
scale extensively with system size, i.e. number of single-particle modes n [15 j. This is 
the reason why in some quantum complexity applications of perturbation theory [2] it is 
required that g\\V\ \ <C A implying an unphysical scaling of g ~ 1/n as ||V|| ~ n. Of course, 
an extensive scaling of errors with computation size, is entirely natural from the point of view 
of a quantum circuit model in which each gate has a certain error rate so that the total error 
rate scales with the total number of gates. Furthermore, in many physical implementations 
of quantum computation, two qubit gate interactions are derived using perturbation theory 
from more basic interactions and even though these qubits all need to couple together and 
thus perturbation theory should be applied at the level of a many-qubit system, this is 
never undertaken. The answer to error accumulation in the quantum circuit model is the 
use of quantum error correction and fault-tolerance so that error rates on logical qubits 
are suppressed [3Tj. We discuss some aspects of making the computation fault-tolerant in 
Section VI. 
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FIG. 1. (a) A pairwise Hamiltonian induces coherent hopping between sites 2 and 3 with amplitude 
g. A pairwise blocking Hamiltonian raises the energy of a particle at site 3 to A when there is also 
a particle at site 1. In degenerate perturbation theory, the sum of these pairwise Hamiltonians 
results in a low-energy sector with an effective three-local Hamiltonian: the presence of a particle 
at site 1 energetically blocks a second particle from hopping from 2 to 3. (b) Quantum field 

effect transistor: the particle at the ‘source’ site 2 can propagate to the ‘drain’ site 4 only if the 
intermediate site 3 is not blocked. Tunneling through the blocked site can be suppressed by adding 
additional intermediate sites with high energy. 

II. HAMILTONIAN QUANTUM COMPUTATION 

Here we describe in detail the execution of a one-dimensional quantum circuit using 
a two-dimensional time-independent Hamiltonian with nearest-neighbor interactions. For 
clarity of exposition, we present our analysis in terms of hopping of spin-1/2 particles. Each 
spin-1/2 particle implements a ‘chronon’ that propagates through the computation, carrying 
with it an internal qubit on its spin degree of freedom. In section [V] we discuss what the 
terms in the Hamiltonian look like when one represents the spin-1/2 particles using dual¬ 
rail qubits, so that the state 100) represents the absence of a particle, 110) represents the 
presence of a particle with internal qubit state 10), and |01) represents the presence of a 
particle with internal qubit state |1) (the state 111) is unused). In the spin-1/2 hopping 
model, particles hop from site to site in the presence of neighboring particles. As they hop, 
their internal qubits can change according to SU( 2) rotations. When the path along which 
hopping takes place is determined by the internal state of a neighboring particle - a spin- 
dependent quantum held effect transistor - then the internal qubit of the hopping particle can 
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undergo a controlled rotation or a CNOT controlled by the internal state of the neighbor. 
Although the goal of this paper is to show how Hamiltonian quantum computation can 
be implemented using pairwise interactions between qubits, the construction given in this 
section using pairwise interactions between spin-1/2 particles also represents a viable path 
to implcmentable Hamiltonian quantum computation, using, e.g., spin-dependent electron 
tunneling [32] , 

For concreteness, we consider a rotated 2D grid depicted in Fig. [2j At each site 
i,j G {0,... ,m} on this lattice a particle with spin 1/2 can reside. There is assumed to 
be only one particle per horizontal line so in total 2m + 1 particles: the dynamics of the 
Hamiltonian preserves this property and will let particles move along the horizontal lines, 
hopping from site to site. At the top and bottom of the lattice, i.e. sites (0, m) and (m, 0), 
the particle is therefore fixed in position and undergoes no dynamics. 

Initially, at the beginning of the computation, the particles are all residing on the left end 
of the grid, one on each site. Using nearest-neighbor interactions on this grid, one can embed 
a one-dimensional quantum circuit with nearest-neighbor interactions (see [2(H [25]). Each 
horizontal line then represents a single qubit wire of the one-dimensional quantum circuit. 
Using SWAP gates such circuit can be used to run an arbitrary quantum computation. 
Instead of using this rectangular grid one can imagine a line of particles hopping forward over 
a grid, —a spatial execution of a one-dimensional quantum circuit-. It was noted however 
in [20, 25] that if the embedded quantum circuit is relatively small compared to m so that 
it can be embedded in the expanding region of the grid where particles are gradually added 
at the boundaries, then the forward motion of the particles in the Hamiltonian computation 
is very efficient. In essence the boundary condition imposed by this grid breaks the time- 
reversal of the computation: the string is more likely to move from the boundary to the bulk 
as the number of bulk strings is much larger. Furthermore, both adiabatic and Hamiltonian 
computation are easy to analyze in this geometry. 

The entire Hamiltonian of the system will be H — H string + g\ h op where H stTing is a 
strong term which enforces particles to stay close together and form a connected string. In 
other words, -String gives a penalty to particle configurations which do not form a connected 
string over the lattice as in Fig. [2j H str ing will have a degenerate ground space with zero 
energy separated by a gap A to higher excited ‘broken string’ states. In H stTing each particle 
interacts via a strong Ising-like interaction with four other particles, namely its neighbor 
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particles on the rotated grid. The perturbative parameter g -C A and 14o P is a sum of 
hopping (kinetic) terms which move particles forward or backwards from site to site (on a 
horizontal line). 

If we treat Vh op perturbatively, the effective Hamiltonian in 1st order perturbation theory 
equals H e 3 = P s trmgI4o P P s trmg + 0(g 2 /A). This effective Hamiltonian is thus comprised 
of hopping terms which preserve the connectedness of the string. Higher-order terms in a 
perturbative expansion correspond to multiple hops of one or more particles which, taken 
together, preserve the connectedness of the string. 

We will first show how one can execute a simple circuit comprised of single-qubit gates 
with H. In order to include gates such as a CNOT (or controlled-P gate), some of the sites 
on the grid will be replaced by pairs of sites, for example, one sitting below the 2D plane 
and another one above. At such locations, the string has the possibility of splitting and 
running through either of the sites, below or above the plane, see an example in Fig. [3] We 
will modify P s trmg so that depending on the internal (spin) state of a control particle, only 
one of these choices has zero energy. Thus P str i ng has zero energy for all connected correct 
strings where the correctness depends on the state of some control particles. Similarly, 14o P 
is modified to allow particles to hop to the newly defined sites above and below the plane and 
undergo internal dynamics which is different for when they hop to a site below or above the 
plane, thereby allowing controlled-rotations and CNOTs. The string is thus being routed in 
two different ways capturing the idea of a railroad switch which temporarily routes a train 
along two different paths. During the quantum computation while the control particles are 
in superpositions of different internal states, one thus works with a coherent superposition 
of strings routed partially above or below the plane. 

The standard procedure to analyze the dynamics of the (space-time) circuit-to-Hamiltonian 
construction is to define a unitary transformation which transforms away the internal dy¬ 
namics of the spin degrees of freedom leaving only the dynamics for the pointer or multiple- 
pointer/string degrees of freedom, see e.g. [T2J. The difference with the previous construc¬ 
tions is the coupling between the internal spin degree of freedom and position of the string 
used in the CNOT or Toffoli gate. In Appendix [A] we show how to rotate away the internal 
dynamics, both for Hamiltonian computation as well as for the quantum adiabatic model 
introduced below. It then follows that the string dynamics due to H e ® is easy to represent 
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as one can parametrize a string in Fig. [2] by a bit-string z of length 2m and Hamming weight 
m, where bit 0 = / (resp. 1 = \). For example, the initial string on the left |z init ) is the 
configuration 100 ... 01... 1). The dynamics of the string in the Hamiltonian computation 
(resp. adiabatic computation), and therefore the forward motion of the computation, can 
then be unitarily mapped onto a one-dimensional XY model (resp. Heisenberg model) and 
this dynamics has been considered previously in [201 [25]. Given that one initializes the 
computation with a string on the left, the probability of finding the string elsewhere on the 
lattice after a time-evolution expis also known [55]. 


The Hamiltonian that we present assumes that there exists only a single-particle or 
chronon on each line. It is not possible to enforce this condition locally, but we will show at 


the end of the paper, Section VIA, that one can add simple local interaction terms (which 
translate as two-qubit ZZ terms in the dual-rail encoding) which ensure that states with 
multiple particles on a line are of higher energy. 



FIG. 2. Sites are labeled (i,j) with i.j € {0,... , m}. The total number of particles is 2 m + 1 
counting the stationary dummy particles at the top and bottom of the grid. There is a single 
particle on each horizontal line of sites. The red string denotes the locations of the particles and 
represents a partially-completed computation. The string can be described by a bit string z of 
length 2m, with Hamming weight m, when we identify / = 0 and \ = 1. 
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A. Single-Qubit Gates 


It is useful to define the particle number operator at site (i,j) as n [i,j] = ^ s=0 1 n s [i,j\ 
where n s [i,j] is the number operator for the particle at site (i,j) in internal state s = 0,1 
(e.g. spin 'l and J.), that is, n s [i,j] = al[i,j]a s [i,j] where a\[i,j] (a s [i,j]) is the creation 
(annihilation) operator for a particle in internal state s at site (i,j). Sometimes we also 
write n[u] or a[v\ etc. where the v label just stands for site v, to avoid cumbersome notation. 
We define 


i/« 


A 


string 


T, h ‘ 


//, 


boundary i 


(3) 


where e labels all the edges e of the 2D grid in Fig. [2j //boundary is a boundary term which 
acts only on the particles at the boundary of the lattice. The goal is for H str ing to have a 
degenerate ground space where each state corresponds to a connected string of particles over 
the lattice as in Fig. [2] It is important to remember that we always work in the subspace 
where there is a single particle per horizontal line. We take an Ising-like interaction per edge 
e: 


H e = -(/ - 2 n [v])(I - 2n[v'}) + /. (4) 

A state with 2 particles at v and v’ (n[u] = n[i/] = 1) or a state with no particles at v and 
v' (n[u] = n[u'] = 0) have energy equal to zero with respect to ^H e and a configuration 
with n[u] 7 ^ n[i/] has energy A. A particle at a site v in the bulk of the lattice participates 
in four edge terms H e . For a connected string going through this site, two out of the four 
terms will have zero energy while the other two edges together give an energy penalty A. 
However at boundary sites, which connect to only 3 other sites, this penalty for a connected 
string becomes A and for the 4 corner sites of the grid a connected string going through 
these sites picks up no penalty at that site. Thus in order for all connected strings to have 
equal energy with respect to H strmg , we need to add A penalties at the boundary sites. One 

takes //boundary 


//boundary = — ( ^ 11 M + n M I ' ( 5 ) 

\ueboundary uecorners J 

A connected string will then have energy (2 m + 1)A. For a broken string at least two sites 
will be end-points of a string: such string thus has energy A above the ground space of 
connected strings. This establishes that the gap of H stl[ngs is A and the degenerate ground 
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space consists of connected strings (independent of the spin-state of the particles) with 
projector P string onto this ground space. 

To implement some single-qubit gates U one takes I4o P = ^hop w h ere V runs over 
plaquettes p of the grid and each plaquette corresponds to some single qubit gate U or I. 
The term Vjf op for a plaquette with corner sites (i,j), (i + 1, j), (i + 1, j + 1), ( i,j + 1) lets 
a particle at site (i,j) hop to site (i + 1, j + 1) while it is changing its internal spin state 
according to the gate U, i.e. 

Ko P = ~ + 1,3 + 1 ]a s [i,j] + herm.conj.. (6) 

s,s'=0,l 

This form of executing single-qubit gates has been first proposed in [22] - When U = I gate, 
this interaction is a simple kinetic energy term of moving particles/fermions on a ID line. 
If we treat H str i ng + gVh op perturbatively, the effective Hamiltonian reads 

H’ s = Jft, ril , i 14 op P„„ tag = g Y. hop + O (d) , 

P ' ' 

^c P ond.hop = - n ^ + 1,3] n[b J + !] ( I^1 «)“!'[* + 1 ’j + l\ a s[hj] + herm.conj. j .(7) 

\s,s'=0,l J 

In words: particles can only hop forward or backward over a plaquette when there are 
particles at the top and bottom of the plaquette so that hopping keeps the string connected. 
It has been shown in [20, 25] that the string dynamics induced by such effective Hamiltonian 
H e ff can be unitarily related to an XY model. 

In principle, it is possible to give an effective Hamiltonian in the string subspace which 
includes higher-order virtual hopping processes from string subspace to high-energy discon¬ 
nected string subspace and then back to string subspace. Such effective Hamiltonian in 
the low-energy string subspace realizes the proper logical single-qubit gates, but the string 
dynamics will not be described by an XY model. For example, in second-order in g, one 
has the following possible processes: ( 1 ) twice application of some V^ op leading to no string 
motion, or a term proportional to /, (2) the application of two different terms Vj^ p and V£ op 
with q 7 ^ p, each of which keeps the string connected leading to terms which are unitarily 
equivalent to Hj- Y , and (3) the application of two adjacent terms hjf op and V^ op with q ^ p, 
such that the application of either one breaks the string but together they map back onto a 
connected string. This third process has a prefactor g 2 /A and does not lead to terms which 
are equivalent to H YY . 
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FIG. 3. The red string denoting the current state of the computation runs from top to bottom of 
the lattice through a set of sites. Two light grey CNOT regions, each comprised of four plaquettes,j 
are depicted. At the center site of these regions the string can either go below or above the plane: 
which path the string is allowed to take depends on the state of the control qubit of the CNOT 
gate. 

B. CNOT Gate 

We describe the modification of the Hamiltonians i/string and gVh op which one allow one 
to perform controlled-C gates. We illustrate the construction with a CNOT gate. In order 
to execute a CNOT we will modify a region of four adjacent plaquettes on the lattice; we call 
this region a CNOT region, see Fig. [4j The central site (i , j ) in this region will have a doubled 
state space; there is a site (say) above and below the plane. The (annihilation) operators 
for this site thus have one additional new label, i.e. they are a s [i,j, k = 0,1] where the k 
label is the new coordinate. In order to give the modified Hamiltonian, we can use number 
operators n a [i,j,k] = k]a s [i,j, k\, n [i,j,k] = E s n s[h^ k l n [bi] = E fc =o,i n [bi, 

etc. We can imagine placing such a CNOT region in various places on the grid. 

We modify the edges which are internal to the CNOT region in H stTing : these are the 
edges a,b,c and d in Fig. |4](b) . We choose the control qubit above the target qubit. The 
control qubit is just undergoing the / gate when hopping from (i, j — 1) to (i + 1, j) in Fig. |dj 
The qubit on the horizontal line below the target qubit and the qubit on the line above the 
target qubit are not participating in the CNOT gate and could be undergoing single-qubit 
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FIG. 4. Left: Grey CNOT regions in the rotated grid. In a CNOT region terms in the Hamiltonian 
are modified in order to execute a CNOT. A long CNOT region is shown at the bottom. Right: 
detailed picture of the basic CNOT region. In the middle of the CNOT region the state space 
is doubled (the target particle can hop to a site below or above the plane). Depending on the 
internal state of the control particle at site ( i,j — 1) and (i + 1, j) an energy penalty is assigned 
when the target particle is present on the wrong side of the plane. Edges which involve these two 
central sites at (i,j, 0) and ( i,j , 1) are labeled a,b,c and d. The horizontal hopping terms labeled 
by plaquette p\ and p 2 towards the doubled [i, j, 0/1) site are modified to execute the proper logic 
on the target qubit. 

dynamics or participate in another CNOT gate. Thus the edges c,d are the same as in 
Eq. (|4j), e.g. H c = —(I — 2n [i — 1 ,j])(I — 2n [i,j]) + I where n [i,j] is total number operator 
of the doubled site (i,j) defined above. 

Instead of Eq. Q, one has for edge a and b 

H a = -(/- 2n s=0 [i, j - 1])(/ - 2n[i,j, k = 0]) 

-(/ - 2n s=1 [i, j - 1])(/ - 2n[i,j, k = 1]) + 21, 

H b =-{I- 272,3=0 [i + 1, j])(I - 2n [i,j, k = 0]) 

-(/ - 2n s=1 [i + l,j])(I - 2n[i, j, k = 1]) + 2J. (8) 

What is the spectrum of, say, H a l Remember we assume that there is at most one particle 
per horizontal line. When there are no particles on sites ( i,j — 1) and ( i,j,k ) for k — 0,1 
H a has zero energy. When there is a particle at ( i,j — 1) but no particle at (■ i,j,k = 0) 
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or ( i,j,k = 1), then jH a equals A/2. When n [i,j,k = 1] = l,n s= i[i,j — 1] = 1, then 
n s= o[i,j — 1] = n [i,j,k = 0] = 0, hence the energy is again zero. When n [i,j,k = 1] = 0 
but n s= i[i,j — 1] = 1, we have n s=0 [i,j — 1] = 0 and two choices namely n [i,j, k = 0] = 0 
or n [i,j, k — 0] = 1. In the first case, the energy is zero, the string simply does not yet run 
over the edge a. In the second case, the energy of ^H a is y, and we can say that the string 
runs incorrectly over edge a: even though the control qubit is in state s — 1, the target 
particle is found at site (■ i,j , k — 0). These arguments imply that not continuing along edge 
a has the same energy penalty as incorrectly continuing along edge a. The same arguments 
apply to the modified edge b. This means that the ground space of H string is formed by all 
correct strings which have energy (2 m + 1)A as before. We denote the projector onto the 
ground space of H stTing as P correc t string- What is the energy of broken or incorrect strings and 
thus the gap of the modified string Hamiltonian? Again a string can be broken such that it 
has A energy above the ground space. A string can also be connected but incorrect on the 
central site in the CNOT region: in that case the string is picking up an energy penalty only 
on the incorrect edge which costs The gap of the modified H stling is thus A/2 (which we 
can of course rescale to A again). 

In order to execute a CNOT gate in a region, the weak hopping terms towards the 
central now doubled site (i,j) also need to be modified. If the particle passes through the 

site (i,j,k = 0), its internal state should be unchanged as the control particle was in the 

state |0). If the particle passes through the site ( i,j , k = 1), its internal state should undergo 
a bit-flip X (for a controlled-!/ gate this can be an arbitrary unitary U). We thus modify 
Eq. (|6]) for plaquettes p\ and P 2 as in Fig. |dj We take 

= — X^s=o i a s[bi) 0]a s [i — 1, j — 1] + herm.conj. I gate when routed through 0 

— i a Mb Ji l] a s[* — 1, j — 1] + herm.conj. X gate when routed through 1 

Ko V = ~ E s =o,i Efe=o,i + 1 O' + 1 \a 8 [i,j, k] + herm.conj. (9) 

With the introduction of the CNOT region(s), the effective Hamiltonian in Eq. (Jt|) , 
H eff = .9 /Correct string 14 op P CO rrect string, is again a sum over controllcd-hopping terms gH? ond hop 
with the modification that for all plaquettes p\ and P 2 in a CNOT region, the controlled- 
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hopping terms equal 


■^cond. hop = ~ n s=o[i,j -l]n[i-l,j] E 0]a s [i — l,j — 1] + herm. conj. 

\s=0,l 

-n s =i\i,j - l]n[i - 1 ,j\ I Y l]a s [i -l,j-l] + herm. conj. 

\s=0,l 

a\[i + 1, j + 1 ]a s [i,j, k] + herm. conj. 1(10) 

s=0,l / 

These terms represent the execution of the CNOT gate via two steps p\ and pp- the target 
particle hops from site (i — 1, j — 1) onto either site (i, j, k = 0,1) depending on whether there 
are particles right above and below it and depending on whether the control particle above 
it is in the state 0 or 1. In the next hopping process P 2 the particle does not undergo internal 
dynamics and it can hop from either site (■ i,j , k = 0) or (■ i,j , k = 1) onto (i+1, j + 1), at least 
when there is a particle in the correct spin state above and a particle in an arbitrary state 
below it. Of course, the Hamiltonian also includes the reverse conditional hopping process. 
One can note that these controlled-hopping processes H^ nd hop and H^ nd hop will take place 
sequentially as the particles above and below the target particle have to move forward or 
backward after one of the process takes place in order to have the next one executed. 

We can consider what terms are present in the next order of the perturbative expansion 
where contributions are coming from twice hopping from the correct string space back to the 
correct string space. If double hopping involves a single particle, then in order for the string 
to stay connected (and correct), the particle has to hop backward and forward, inducing no 
effective dynamics. The only other contributions in second-order in the perturbation Wop, 
are from two different particles hopping so that together they preserve the connectedness of 
the string. The same is true for all third-order processes as it takes at least four hopping 
processes in order for the particles to hop/tunnel through a CNOT region. Thus terms up to 
0(g 3 / A 2 ) give rise to a change in string dynamics, but the logic is still properly executed. In 
other words, if the string has traveled through the region where the gates have been applied, 
then the correct circuit has been implemented. 

A fourth-order process can however lead to incorrect logic in the CNOT gate as follows. 
Assume the control particle is in the state |0). The target particle hops from site (i — 1, j — 1) 
through incorrect site (■ i,j , 1), and then again to (i + 1, j + 1) and the control particle hops 
from (i,j — 1) to (i + 1, j), while the particle below hops from (i — 1, j) to (i, J + 1), so that the 


TTP2 

cond.hop 


n k [i + 1, j]n [i,j + 1] 


fc=0,l 
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string stays connected and is correct after these four processes. In this fourth-order process 
(of order 0(g A / A 3 )) the target particle undergoes a bitflip although the control qubit was 
0. One can suppress these higher-order virtual tunneling processes by replacing the short 
CNOT by a long CNOT. The idea of the long CNOT is to enlarge the CNOT region to an 
elongated CNOT region in which the central site is a sequence of L sites where the string 
can only go below or above the plane. The hopping terms of these L sites propagate the 
particles on one site of the plane or the other. For such long CNOT with L central sites, it 
takes L + 1 hops for the target particle to be on the plane again where the fact that it has 
taken the wrong route is no longer visible (in other words, it takes L +1 hops for the particle 
to tunnel through the energy barrier). For the string to stay connected, the particles above 
and below the target particle also have to hop along, over L sites. Hence, the first term 
which gives rise to incorrect logic is of order O (g(J() 3L+1 ). 

III. TOFFOLI GATE AND PROGRAMMABLE CIRCUITS 

The principle of the spin-dependent quantum held effect transistor can also be used to do 
a Toffoli gate or make the quantum circuit programmable. In the Toffoli gate two particles 
determine how a third qubit should be routed, above or below plane. One can choose these 
two particles as a particle above and below the target particle, as in the CNOT region 
in Fig. [4j For a short Toffoli gate, the modification of H stimg then involves all four edges 
a,b,c,d, each one of them should depend on the internal state of the control particle(s), as 
in Eq. (|8]). 

In order to make the computation programmable, one can add classical control bits to the 
grid which can route the particles to regions where they can undergo the desired logic. These 
bits do not need to participate in the dynamics, i.e. but they can interact with the particles 
involved in the computation in the Hamiltonian H string giving penalties for incorrect string 
configurations. 


IV. QUANTUM ADIABATIC COMPUTATION 

To execute an adiabatic quantum computation on the grid, we define the Hamiltonian 
H (A) with A € [0,1]. This construction follows the ideas in [25], but the adiabatic Hamil- 


19 


tonian has to be chosen differently in order for the spectral analysis of [25] to apply. Fur¬ 
thermore, for adiabatic computation the previously defined Hamiltonian H = String + 9 Wop 
does not suffice: new interactions have to be included in the adiabatic Hamiltonian and the 
adiabatic parameter has to be chosen so that the ground state is a (weighted) superposition 
over all strings, thus encoding the computation. 

We will take //(A) = -Hiring + //input + g //circuit(A). Here //input is a term setting the 
initial internal state of the particles to the correct input state of the circuit: one can simply 
choose the proper occupation of sites on the left half of the lattice. If these initial states 
are set correctly, // input has zero energy and thus keeps the degeneracy between all correct 
strings. We will further neglect //input in the discussion here (see j25]J for how //input makes 
the ground state of the quantum adiabatic computation unique). 

We construct //circuit(A) as a modification of Wop- Let us first do this for all plaquettes 
which correspond to single-qubit gates. Let a plaquette p enable a particle to hop from site 
(i, j) to ( i + 1, j + 1) while there is a particle at the top of the plaquette, site (i + 1 ,j) and 
at the bottom of the plaquette (i, j + 1). The idea of the quantum adiabatic computation is 
to gradually turn on the motion/kinetic energy of the string using a parameter A. Initially 
for A = 0, we want to choose //circuit(A = 0) such that the straightest string, with the fewest 
‘wiggles’, has the lowest energy. There are two such straightest strings on the left and right 
boundary of the lattice and we include a term H init in H circuit (A) which favors the initial 
string on the left. We choose H in it = n[m, 1] + n[l, m], thus penalizing any string which runs 
along the top and bottom edges on the right. 

In the string subspace, it is clear that a term of the form n [i + l,j]n[i,j +1] provides an 
energy penalty when the string runs around the plaquette p of which (* + 1 ,j) and (■ i,j + 1) 
are the top and bottom sites. More precisely, we will take 


//circuit (A) = ^ ( n t* + lj /] n [b / + X ] + XV hop) + y/l- A 2 H init , (11) 

P 

where is as before. The additional terms that we have thus added are diagonal in the 
string basis. Thus in the string subspace, the effective Hamiltonian just equals 


H&mW = 9 ( Y1 t j]n[t j + 1] + + Vl - A 2 // mit , 


( 12 ) 


When the quantum circuit contains CNOT gates, we can define the adiabatic Hamiltonian 
in the same way, using Eq. (11). For certain plaquettes p in the CNOT region V^, is chosen 
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according to Eq. ([9]) and for a doubled-site (i,j) we simply use the definition n [i,j] = 
n [i,j, k], that is, we use the total particle operator on the site. 

It can be shown for this effective Hamiltonian H^ CVLit ( A) that one can rotate away the 
internal spin dynamics including the dependence on the position on the control qubits, see 
Appendix |A} It then follows that we can represent the action of the effective Hamiltonian 
by a XXZ Heisenberg chain as in [25] in the string subspace spanned by 2m-bitstrings | z): 

- 2m—1 

( Z |^(n[i + 1 ,i]n[i,i + 1 ] + A< nd . hop )|z> = -~(z\ (Z j Z t+1 -/)-A(A,A' j+1 + y i y i+1 )|z}. 

p i= 1 

(13) 

In addition, the term H m - ]t can be represented as a kink boundary condition which allows 
this model to be exactly solvable: thus all previous results concerning the ground state and 
the gap in [25] apply to this effective Hamiltonian. 

V. DUAL RAIL ENCODING: REQUIRED INTERACTIONS 

There are different ways of representing the particle with spin physically: the simplest is 
a dual-rail encoding using two qubits to represent one spin-1/2 particle. We note that the 
fermionic nature of the particles in the construction is not relevant: each fermionic mode 
can be represented by a qubit where occupation of the mode is |1) and no occupation of 
the mode is |0) (see [33] for the argument). This means that the two fermionic modes of 
a particle with spin can be represented by two qubits. We label these two qubits as qubit 
0 and 1. We can choose to represent a particle in state s = 0 as |10)oi, he. qubit 0 is 11), 
particle in state s = 1 as |01) 0 i, no particle present is |00) O i and the two-qubit state |ll) 0 i 
is not used as there is at most one particle present. 

Thus, each creation operator a\[i, j, k] (or annihilation operator a[i, j, k ]) with spin-label 
s can be represented as a single-qubit operator cr+[i, j, k] = |l)(0| s [i, j, k] (resp. k] = 

|0)(l| s [z, j, k]). This implies that all hopping terms in I4 op , Eq. (§. are weak two-qubit 
terms of the general form 

I or X gate - (cr+crj" + cr+) = — ],(A',Xj + YjYj) 

Hadamard gate |^=( X. L Xj + k/Y)) 

complex single qubit U gate oc X/Xj + Y t Yj and oc X'iYj + XjYj 

Toffoli or CNOT gate - ^(X l X J + Y^) (14) 
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In Fig. gc ) we show these interactions for some specific single-qubit gates. 

What about the diagonal terms in Z/ string (or H cilcuit (X), Eq. (11))? 

One has Z 0 \v\ = I — 2n s= o[v\ and Z x [u] = / — 2n s= i[v\ for any site v. Using this identih- 
cation, a normal edge H e in //string as in Eq. Q then becomes the Ising interaction 


He = ~{Z Q [v] + Z^v] - I)(z 0 [v'] + Z^v'] -/) + /, 


(15) 


hence four pair-wise Ising interactions between the four qubits on sites v and v' (and a sum 
of local Z terms on all 4 qubits). The boundary terms acting only on qubit sites at the 
boundary of the lattice, Eq. ([5]) can also be represented as single-qubit Zq\v\ and Z\[v\. A 
modified edge term in a CNOT region, say, H a in Eq. (J8]), simply equals 

H a = ^Z Q [i,j-l](Z o [i,j,0] + Z 1 [i, j,0]-/)-Zi[i, j-l](Z 0 [i,j, 1] + Z 1 [i,j, 1}-I) + 2I. (16) 

For adiabatic quantum computing the Hamiltonian in Eq. ([TTJ) contains additional inter¬ 
actions which can be represented as weak Z and two-qubit ZZ interactions. In Fig. [5] we 
show the connectivity of the two-qubit interactions for Hamiltonian and adiabatic comput¬ 
ing. One can note that qubits connected via Ising ZZ terms are never connected via XX+YY 
or XY+YX terms and vice versa. For Hamiltonian computing, each qubit 0 or 1 at a site v 
in the bulk of the lattice interacts with 8 other qubits, along the four different edges, via a 
ZZ interaction. When that qubit is part of the control particle of a CNOT gate, it interacts 
with 10 other qubits as the state-space is doubled on the central CNOT site. 

If we assume that we intersperse single-qubit gates and CNOT gates with I gates, then 
each qubit interacts in addition with at most 3 other qubits via XX+YY and XY+YX 
interactions. For adiabatic computing each qubit at a site interacts in addition with 4 other 
qubits (above and below the site), via ZZ, see Fig. [5](b). 

A possible universal gate set is comprised of the T = diag(l, e 77 ? 4 ) gate, a CNOT gate and 
a Hadamard gate. For the T gate one needs to connect two qubits via a linear combination 
of XX+YY and XY+YX, see Fig. gc ), which requires more physical interaction engineer¬ 
ing than only XX+YY. Alternatively, one can use Hadamard and the Toffoli gate to get 
universality [ : 34j. In this case both for the Hamiltonian computation and the adiabatic com¬ 
putation, one only needs the interaction ZZ, Z and ± XX+YY. The set of interactions ZZ, 
Z and -(XX+YY) are stoquastic (sign-free) and can be mapped via perturbative reductions 
onto a transverse field Ising model [29j. Note that in order to do the Hadamard gate one 
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needs ± XX+YY and thus the quantum adiabatic Hamiltonian which realizes a universal 
computation using a Hadamard and Toffoli gate is not stoquastic or sign-free and is unlikely 
to map onto a transverse held Ising model quantum annealer. 




(a) (b) extra for adiabatic QC 

— ZZ — XX+YY — XY+YX 



I gate X gate Hadamard S gate T gate 


(c) 

FIG. 5. Required two-qubit interactions for Hamiltonian and adiabatic quantum computing. At 
each site the particle with spin is represented by two qubits (0 and 1) in the dual-rail encoding, 
(a) For each edge connecting a pair of sites there are 4 strong ZZ interactions (in red) between the 
four qubits. Thus each qubit interacts with 8 other qubits via ZZ. (b) Only for adiabatic quantum 
computing there are, in addition, weak ZZ interactions between qubits: for each plaquette the two 
qubits at the top interact via ZZ with the qubits at the bottom, (c) Single-qubit gates are realized 
by a plaquette and require weak interactions XX+YY and/or XY+YX for qubits on the left and 
right of the plaquette. Here S = diag(l,i) and T = diag(l, exp(i7r/4)). For each plaquette in the 
total grid which is not used for a non-trivial gate, one has an I gate. The CNOT or Toffoli gate 
(not shown) use sites above and below the grid, similar strong ZZ edges and an X and I gate. 


A. Physical Implementations 

The methods for Hamiltonian quantum computation developed in this paper are po¬ 
tentially suitable for implementation on multiple physical platforms. In the form given 
here, we require (a) strong ferromagnetic ZZ interactions between neighboring qubits (plus 
strong local Z terms), and (b) weak XX+YY interactions. For example, superconducting 
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transmon qubits may be suitable to realize the 2-qubit interactions for Hamiltonian com¬ 
putation described in this proposal. The strong ZZ interactions can be obtained by strong 
capacitively-couplcd ‘octmon’ (in analogy with Xmon qubits [35]) qubits which have eight 
arms through which they couple with nearest-neighbors octmons. The weak XX+YY in¬ 
teraction between two qubits across a plaquette can be obtained by placing a bus-resonator 
on the plaquette through which both qubits couple. This resonator-mediated interaction 
through which the qubits can virtually exchange a photon is the basis of the iSWAP gate 
[3BJ. The coupling strength J of such term J(af crj + crjcrf) scales as J = fi'jc/A where A is 
the detuning between resonator frequency and qubit frequency (assuming that both qubits 
are at the same frequency) and gjc is the Jaynes-Cummings coupling. The frequency of 
each bus-resonator and thus A can be adjusted to fix the strength of the weak coupling J. 
In such a physical set-up one would realize Hamiltonian computation in a rotating frame of 
each octmon qubit, where the rotating frame is chosen such that one retains single-qubit Z 


terms of the correct strength, as in Eq. (15), (16). 

Adiabatic computation is less suitable for qubits whose spectrum is nondegenerate since 
the thermal environment tends to relax the qubits to their individual ground state |0) instead 
of the joint multiple-qubit ground state of the adiabatic Hamiltonian. Thus other qubits 
with a degenerate spectrum such as flux-qubits currently used in quantum annealing could 
be considered for an implementation of adiabatic computation. 

A direct realization using electrons and their spins may also be possible: at least one can 
imagine a spin-dependent hopping term realizing the Hadamard gate. In this direct realiza¬ 
tion the Hamiltonian requires a strong attractive interaction between neighboring electrons 
propagating over one-dimensional lines (in the adiabatic model there is in addition a repul¬ 
sive interaction between next-nearest neighbor electrons across a plaquette, see Fig. [5]). One 
could consider using an particle/hole encoding alternating for adjacent wires so that parti¬ 
cle and hole on nearest-neighbor wires attract via the Coulomb interaction (while repelling 
each other on next-nearest neighbor wires). Another idea is to use spin qubits localized in 
arrays of quantum dots and obtain the XX+YY coupling by letting two spin qubits virtu¬ 
ally exchange a boson via a quantum bus. Such a quantum bus can be a superconducting 
microwave resonator as in e.g. [37] or a standing surface acoustic wave as proposed in e.g. 
J. Strong ZZ coupling could again be obtained via an electrostatic interaction. 

We emphasize that these proposed implementations are simply sketches of how to attain 
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the requisite interactions using existing technologies. Any quantum technology that allows 
the implementation of strong Ising couplings and weak hopping terms allows the construction 
of a Hamiltonian quantum computer in principle. In the next section we discuss various 
obstacles that stand in the way of constructing such quantum computers in practice. 


VI. INACCURACY, NOISE AND ERRORS 

Any physical implementation of the construction given here will be susceptible to man¬ 
ufacturing inaccuracies and to dynamic noise, which will induce errors in the computation. 
Instead of simulating a one-dimensional quantum circuit the Hamiltonian or adiabatic com¬ 
putation could simulate a fault-tolerant quantum circuit. In general such fault-tolerant 
quantum circuit, for example the surface code architecture EB, is a 2D nearest-neighbor 
circuit. This implies that the time-independent Hamiltonian or adiabatic computation takes 
place on a 3D grid where the degrees of freedom are represented by a 2D-dimensional prop¬ 
agating membrane. Quantum measurements of subsets of qubits during the computation 
(to implement quantum error correction) can in principle be incorporated in any geometry 
by terminating the spatial regions of qubits which are to be measured so that they cannot 
propagate any further. No rigorous mathematical analysis of adiabatic or Hamiltonian com¬ 
putation on such a 3D grid with regions in which partial membranes terminate early has 
been undertaken, although one expects similar (spectral) results to hold as in the 2D grid 
case. 

Running the encoded quantum circuit could provide protection against all errors or in¬ 
accuracies in the Hamiltonian or adiabatic computation which can represented as errors on 
the qubits in the simulated fault-tolerant circuit as long as the induced errors have sufficient 
locality in space and time |39|. We discuss a few issues to be resolved. Consider what 
happens if the string dynamics is very slow so that a single-qubit noisy gate realized in some 
plaquette is repeatedly done and undone. It means that the error rate on this single-qubit 
gate will be high. It was argued in [25] that the forward motion of the string is constant 
in the expanding region of the grid, but does this still hold for a noisy version of the string 
motion? 

Static disorder in the Feynman Hamiltonian leads to Anderson localization of the clock 
variable [JO]. While the extended nature of the multi-clock string may render the Margolus 
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Hamiltonian less susceptible to localization, no analysis of localization in this system has 
been performed. 

Static or dynamic disorder in the ideal time-independent Hamiltonian or coupling to a 
finite temperature bath also leads to errors which represent leakage errors in the original 
encoded circuit. For example, in the dual rail encoding, amplitude damping on one of the 
qubits at a site leads to the state 100), i.e. the loss of a particle. Bitflip errors on the pair 
of qubits at a site can lead to the state 111) and thus again a leakage error. 

Another important source of errors (due to perturbative corrections, static disorder in the 
blocking terms or thermal activation of forbidden tunneling) leads to the string to become 
disconnected or incorrect. This can lead to the computation getting stuck: the forward 
motion of a particle 1 is blocked as the particle, say, 2, which gates this motion has already 
moved further in its computation. However the forward motion of particle 2 will continue 
to be blocked while its backward motion, rejoining particle 1, can lead to a lowering of 
its energy. Incorrect string states in which by some means a particle gets on the wrong 
CNOT track and picks up a X can be surpressed by making a CNOT gate long. One can 
thus expect that a sufficiently low-density of string breaks or wrong string turns can be 
dealt with by simulating a fault-tolerant quantum circuit which includes leakage protection. 
Leakage errors are of course also an issue in circuit-based quantum computation requiring 
leakage reduction units such as quantum teleportation or swapping into fresh qubits (see 
e.g. [?TJ and references in [[31]). 

A. String Loops and Multiple-Time Wavefronts 

Last but not least, we consider what happens when the number of particles on a line 
is not conserved, or in the dual-rail encoding, bit-flip errors generate new particles on the 
line. These are dangerous sources of errors which in the presented model can lead to lower 
energy states if we do not add additional penalty terms. If there is at most one particle per 
line the string always runs down the lattice as in Fig. [2] However, if one allows for errors 
which create new particles, then one can create a string which loops back onto itself. The 
terms of // s trmg are such that when one particle is present on a site, it is favorable to have 
all nearest-neighbor sites be occupied with particles, so a string which comes from the top, 
loops back and closes onto itself will in fact have lower energy than the ground space of 
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^string in the single-particle per line sector (as on the looping string there is one site where 
a particle has 3 neighbor particles, elsewhere each particle has 2 neighbors). But we can 
add a term H no i oop = E ]Tb . n[i, j]n[i + 1, j + 1] to the Hamiltonian where E S> A which 
gives an energy penalty for a state which has two particles next to each other on a line. 
This term gives a direct energy penalty to a string looping back. The term does not prevent 
there from being another piece of string, i.e. a set of particles representing a new partial 
wavefront, that is created further away on the grid. However, an open string segment will 
pick up energy penalties at its boundaries due to H stTing and a closed string segment will 
pick up energy with respect to H no i oop . 

Energetically penalizing such errors is particularly relevant in the dual-rail encoding where 
the particles are represented by pairs of qubits and thus bitflip errors can create particles. 
Choosing |00) as the lowest-energy no-particle state is useful so qubit relaxation will not 
cause the creation of new particles. For qubits, the term H no loop translates into a set of 
ZZ terms between the qubits horizontally across plaquettes. ft is an open question whether 
the physics that one obtains when such errors are suppressed but not eliminated, e.g. the 
splitting, creation and collisions of additional (closed or open) wave fronts, can be handled 
by error correction. 

One expects that the membrane computation in 3D will be more robust than string 
computation in 2D. A piece of open membrane which has been disconnected from the rest 
of the membrane computation costs an energy scaling with the perimeter of the membrane, 
thus suppressing the formation of large open disconnected membranes. Closed membranes 
are made energetically unfavorable by using no-loop terms as in the string computational 
model presented in this paper. 


VII. CONCLUSION 

This paper has shown how to perform universal Hamiltonian and adiabatic quantum com¬ 
putation using only two-qubit interactions and without higher-order perturbation gadgets. 
Instead, energy penalties are used to enforce the desired Hamiltonian at lowest order. The 
resulting low-energy Hamiltonian could be used either in a dynamic fashion or in an adia¬ 
batic/quantum annealing context to encode the quantum computation in the ground state. 
The construction given here was designed to be potentially physically realizable using large- 
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scale quantum integrated circuits consisting of superconducting qubits or quantum dots. 
Many open questions remain, notably how to deal with errors induced by manufacturing 
inaccuracy and by dynamic noise. 
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Appendix A: Rotating Away The Internal Dynamics 


We assume a Hamiltonian computation on the rotated grid with effective Hamiltonian 
H eff = g Ep ^cond.hop with #cond.hop as in Eq. ([0 inside CNOT regions and Eq. 0 for 
single-qubit gates or / wires, acting in the correct string subspace. Strings in this subspace 
can be parametrized as |z, s) = |z, si = ±1,... = ±1) where 2 is a 2m-bit string 

representing the string degree of freedom and s is a set of k(z) binary labels, one for each 
short CNOT region through which the string z goes, which indicate whether the correct 
string goes below or above the plane. The values s l = ±1 thus correspond to the internal 
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states of the control particles through which the string z goes. Clearly, the number of these 
additional labels k(z) depends on the string z, i.e. some strings do not have any additional 
labels. 

We can write P C orrect string = Ylz,a k s )(^ s l- 

Let V( 2 , s| 2 ', s') be the transformation on the internal states which is executed by the 
circuit in going from a correct string ( 2 ', s') to another correct string ( 2 , s). In other words, 
V is the total matrix-valued amplitude of the path from string ( 2 ', s') to ( 2 , s). We can write 


H eff = - 9^2 f / (^,s|2 , ,s , )|2,s)(2 , ,s'|. (Al) 

(jz,s) (z',s')sNeigh(z,s) 

where Neigh (2, s) are the correct strings which are neighbors of (2, s) meaning that condi¬ 
tional hopping terms in H e & can map (2, s) onto (2', s') and vice versa. 

In addition, we have the composition V(z, s 1 2", s") = ^(2, s 1 2', s')V r (2 / , s' 1 2", s") 

where V(z, s|2", s") is independent of the intermediate string 2': all paths from (2", s") to 
(2, s) execute the same computation independent through which intermediate string 2' they 
go. A transformation V(2|2",s") is unitary on the internal state space if the string 2 has no 
control particle labels, hence we do not post-select the dynamics on the final state of some 
control particles. It also holds that V^z, s|2 7 , s') = V(z', s'| 2 , s). 

We write down an isometry W such that H e ^W is a transition matrix in the string 
subspace \z) independent of the internal dynamics and the spin-labels s: 


W = ^I/( 2 ,s| 2 imt )| 2 ,s)( 2 |, (A2) 

Z,S 

so that W^W = W 1 2)(2 1 = -Paring- Note that W is an isometry as it maps from the string 
space to a higher-dimensional correct string space where the string has additional labels. 
Using that ^) s V(z, s1 2 ', s')V(z', s' |2i n it) = V(z, s 1 2i n it), V(z| 2 ) = / etc., one obtains 




(A 3 ) 


W ] H eff W = -g 

z,z r GNeigh(2:) 

representing only the string dynamics on the grid. 

In the quantum adiabatic computation the effective Hamiltonian, Eq. 0 is slightly 


different than the Hamiltonian in Eq. (Al). It includes a diagonal term of the general form 
s a A z i s )(u s | where the weights a z do not depend on the internal states of the control 
particles s. Conjugation with W of such term gives a z \z, s)(z, s|VU = Yh z a A z ){ z \ 1 
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again removing the internal control particle labels. In addition, the effective Hamiltonian 
includes a diagonal term H input which penalizes the internal states of particles before they 
run through the computation: this means that the strings (z, s) for which this term is non¬ 
zero are always strings z without control-particle labels and V(z, s|zi ni t) = V(z|zmit) = /, 
or the dynamics from z init , to these strings z is trivial. This implies that conjugation by 
W leads to a term which penalizes the internal state of particles before the computation is 
executed, identical to the effect of such term in [25] . 

In all, what these arguments show is that even for a circuit with CNOT and Toffoli 
gates, one can transform away the unitary dynamics of the computation and obtain a model 
which purely represents the string dynamics. One can then invoke the previous results 
which analyze the string dynamics in the case of Hamiltonian computation (XY model) or 
quantum adiabatic computation (ferromagnetic XXZ chain with kink boundary conditions, 


Eq. ©). 

It can be noted that the inclusion of higher-order terms in the perturbative expansion 
(of 0(g k+1 / A k )) in the correct string subspace do not affect the validity of this analysis. 


Such higher-order terms, due to multiple hopping, lead to new terms in Eq. (Al) of the form 
r )z,z , V(Zi s|z', s')|z, s)(z 7 , s' | for some coefficients ~{ zz j where (z, s) and (z 7 , s') are not direct 
neighbors, but are connected through multiple hops. 
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